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I. INTRODUCTORY REMARKS 

A mapping of the spin-i operators onto Fermi operators by means of the Jordan- Wigner transformation was 
used by Lieb, Schultz and MattisS to introduce the exactly solvable one-dimensional spin-i XY model. Later 
the famous Onsager's solution for the two-dimensional Ising model was reproduced by using the Jordan- Wigner 
transformation for the transfer matrix of that modeli. The Jordan- Wigner transformation is the essential ingredient 
of the studies of the statistical mechanics properties of quantum spin chainsi. Much effort has been devoted to 
generalize the fermionization procedure for twdHil and threeSEl dimensions. In the present paper we review the 
Jordan- Wigner transformation in two dimensions as well as some existing applications of this mapping for the spin 
system theory. 

We shall consider a spin model consisting of N — N x N y (N x — + oo, N y ~ > oo) spins i on a square lattice of the 
size L x Ly (L x 

— > oo, L y — > oo) governed by the Heisenberg Hamiltonian 

H = M ■ Sj + hs l 

(iJ) 1 

OO OO oo oo 

i=0 j=0 i=0 j=0 

Here (i, j) denotes all different nearest neighbouring sites at the square lattice, J is the exchange interaction between 
the neichhonrinp' sites, h is the external field. The isntronic Heisenherc interaction consists of the isotronic XY 



part and the Ising part 



S; • Sj 



sf + s»s» +sf S ? 
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(2) 

i. It is important 
(3) 



(4) 

The aim of the Jordan- Wigner trick is to transform the spin variables into pure fermion variables. 

We start the paper by giving a short reminder of the Jordan- Wigner transformation in one dimension (Section 
II). Then we discuss the extensions for two dimensions suggested by M. Azzouzi (Section III), Y. R. Wang! (Section 
IV), and E. FradkinS (Section V). The mean-field-like treatment of the Jordan-Wigner fermions for the isotropic 
XY model is discussed in detail in Section VI. The consideration of a model with Ising term in fermionic language 
is given separately in Section VII. Finally, we summarize some of the results obtained using this approximate 
approach and comment on a comparison with the results derived using other methods (Section VIII). 



where we have introduced the spin raising and lowering operators = s x ± \s y and s z = s + s 
to note that the operators s + , s~ obey the Fermi commutation rules at the same site 

{ S r, s +} = 1, { a +, s+} = { S 7,s-}=0 
and the Bose commutation rules at different sites i 7^ j 

1 = 0. 



II. THE JORDAN-WIGNER TRANSFORMATION IN ONE DIMENSION 

With the help of the Jordan-Wigner transformation we introduce instead of the operators s + , s~ satisfying (^), 
the operators c + , c satisfying the Fermi commutation rules (both at the same and different sites) in terms of 
which the Hamiltonian of the isotropic XY chain is a bilinear form and the Ising interaction yields the products 
of four Fermi operators. Explicitly the Jordan-Wigner transformation reads 

r + — pia« s + c — -i"n„- „+ _ „-ia n + - _ ia„ 

n-1 

a n =TT^n 3 , nj^c+Cj. (5) 
j=o 

The signs in the exponents and the order of the multipliers (first line in ([|)) are not important. To get the operator 
c n (c^) we must multiply s~ (s+) by the exponent containing a sum of rij at all previous sites < j < n — 1 as 
can be seen from (^|) and is shown symbolically in Fig. 1. 
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FIG. 1. Towards the Jordan-Wigner transformation in one dimension. 
To demonstrate how the transformation (|^) works we need the following equations for Fermi operators 

cje j 3 = cj , e 3 C J = ~ C J ■ c j c 3 = — c j> c 3 c j = c j' 



(6) 



which can be easily obtained since 



eW c , - 1 + mcfcj + i (Itt) 2 cfcj + . . . = 1 - cj Cj + e^cj Cj = 1 - 2cj Cj . 
From (^) one finds that 

e [ *< c - Cl = c ; e i,rc " c ", e i7rc - c "c+ = c+e i7rc " c ™, (7) 

if n 7^ i, but 
Besides, 

gSiTTc+c, = L (g) 

Let us show that the commutation rules for the spin raising and lowering operators s + , s~ (^) are given by (^), 
(^) if c + , c are Fermi operators. At the same site we find 

S n S n = c " e " e ™ c n = c »i c ri I s n s n = c n c nj s n s n = c n c n ~ 0, S n S n — C n C n = 0, (10) 

and as a result Eq. (||) becomes evident. At different sites n and n + m (without any loss of generality m > 0) we 
have 



— 4- i7r c i c i + + i 71 " >, c z Cz 

n S rx+m — C " e ' = ™ C n+m — C « C n+m e ,= ™ I 



s+ s~ - r + P l7T l^i- n c i c, r — -r + r P 1,r L,-„ c i c ' (11) 
A n+m ft )i — Si+m L "--"-»> L n — O n+m O n e "i-n , ^llj 

etc. which immediately yields Eq. (||). 

Let us write down the transformed Hamiltonian. Besides s^sj — c^Cj (|l^) we have 

sjs j+1 =cje i ] c j+1 = cjc j+1 , sj sj +1 = Cje ' 3 cj +1 = -Cjcf +1 , (12) 

and as a result the one-dimensional spin-i Heisenberg Hamiltonian (|l]), (^|) becomes 

F = EQ J (cpi+i - + J (ct Cj - 5) ('V:<V - 5)) + £ h I) ■ ( 13 ) 

It is clear now to what extent the Jordan- Wigner transformation simplifies further statistical mechanics cal- 
culations. A nontrivial in spin language isotropic XY chain transforms into tight-binding spinless fermions and 
further rigorous treatment becomes possible. For the anisotropic XY chain the operators CjCj +1 , CjCj+i enter 
Eq. ( |T~3| ) besides, and therefore the Bogolyubov transformation is required in addition. The Ising term leads to an 
interaction between the Jordan- Wigner fermions, however, the low-energy properties may be analysed using the 
bosonization techniques!. While interesting only in the low-energy physics of s > | spin chains one may represent 
the spin-s operators as a sum of 2s spin-^ operators and then, making use of the Jordan- Wigner representation 
for the latter operators, proceed in fermionic language. 



III. THE JORDAN— WIGNER TRANSFORMATION IN TWO DIMENSIONS (M. AZZOUZ, 1993) 



Consider the spin model (|l|) on a square lattice (Fig. 2). 
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FIG. 2. Towards the Jordan-Wigner transformation in two dimensions! 



Two coordinates i and j are taken at the x and y axes, respectively, to specify a given site. 
M. Azzouz definedi the extended Jordan-Wigner transformation as 



s i,j ~ Ci d e 



la i}j _ la; 



i—1 oc j— 1 



*h3 



= ?r EE^z+E^/ ' '•,/./«./ (14) 

\d=0/=0 /=0 / 

(compare with Eqs. (^)). The signs in the exponents in ( Jl4| ) (and the order of the multipliers in the first line in 
(pT[)) arc not important. Let us show that the introduced transformation (14) enables one to construct a fermion 
representation for two-dimensional spin-i models. 

We start with the commutation rules. At the same site one has 



S q,p S ~q,p ~ C 1'P" 



—la + 

Cn r>e ' e 1 C n n Cn TfC n n , S n n S f< Tl C n ^Cn T?i ^n n n ^n « C n ^, S n y.S^ ~ Cn T)C, 



./-/. ~ ~q,p c q,pi s q,p s q,p ~ c q,p^q,p^ 



+ „+ _ „+ „+ 

q,p°q,p ~ ^q^q-.p-i s q,p°q,p ~ ^QtP^tP) 



(15) 



and the Fermi type commutation rules remain unchanged. To illustrate how transformation ( [T4| ) works in the case 
of different sites let us consider, for example, two sites q,p and q,p + m, m > 0. Then, similarly to Eq. ( p~T| ) one 
finds 



s - s + - c e ^f= P nqJ »+ 

b q,p°q,p+m ~ c q,P K 



C q,p+m — c q,P c q,p+m e f=P ' > 



+ - + 1 

q,p-\-rn^q,p 9-P+' m 



-g,p — °<},p+m <?,P c (16) 

(we have used an analogue of Eqs. ([?]), (||) in two dimensions) and hence the operators s~ p and s^ p+rn commute if 
the operators c 9iP and c^ p+m anticommute. Following the same reasoning one can check the rest of the commutation 
rules. 

Consider further the transformed spin Hamiltonian. Let us treat somewhat more general nearest neighbour 
interactions than the one in Eq. ([!]). Namely, we assume different values of interaction in different directions on a 

srmarp latti rp as shown in Ficr 3. 
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FIG. 3. Nearest neighbour interactions on a square lattice. 

(Evidently, we can perform the fermionization presented below for a completely nonuniform model characterized 
by a set of intersite interactions {. . . , Jij-.i+ij, Jij-.ij+i, ■ ■ ■}■) If J = J' = J± = J'^ one faces the uniform 
square lattice. To examine the effects of interchain interactions in quasi-one-dimensional systems one may consider 
the case J±, J' ± <C J, J' (with J not equal to J' if the dimerised chain is considered). If J' ± <C J, J', J± one has 
a model of interacting two-leg ladders. In the limiting case of J' ± = (the noninteracting two-leg ladders) the 
model may be reduced to a one-dimensional system with interactions extending over the nearest sites. 

We begin with the isotropic XY interaction (see Fig. 4). 




FIG. 4. Towards the fermionization of Hxy- 



Inserting ([u]) into (0), (||) one finds 

_^ oo oo 

HXY — ^ X/ X/ ( 8 i,j 8 t+lj + S ij' S i+l,j) + { S i,j S t,j+l + S iJ S i,i+l)) 

i=0 j=0 



2 

i=0 j=0 



i=0 3=0 ^ ^ 

+ + 1 (-Cijc+- + 1 + C+-C 4J + i)) (17) 

(to get the last equality we have used an analogue of Eq. (||) in two dimensions). After introducing the notations 

(oo 3-1 
E • E "' • 1 * 
f=j f=o 

( oo j-1 \ 

0i,i+lO') = 7T E + E ni+1 '/ ' 

\/=J + l /=0 / 

Pi,i+i(»)='"»i,j (18) 

the Hamiltonian (|l7]) becomes as follows 

oo oo 

H ** = ! EE (cue-^y + ^M^ 

a.j. ..... ( P . . e -ivj.j+i(*) r + i r + . , , V 

. oo oo 
i=0 3=0 

+^i,i;i,3+l (-Ci,jcf j+1 + C+-C 4J + i)) . (19) 

Eq. ([l9]) can be viewed as the Hamiltonian of a two-dimensional tight-binding-like spinless fermions with the 
hopping amplitudes 

T (20) 



2 

i=o 3=0 



in the x direction and 

in the y direction. Those hoppings depend in a complicated way on a configuration of the 'intermediate' sites. Their 
complexity explains how the isotropic XY model becomes difficult to examine in two dimensions in comparison 
with an obvious analysis in one dimension. 

There are no difficulties in rewriting the Ising interaction in fermionic language 

1 



8 i,j+l a i,j+l 2 



EE ( J W+1.J ('•-'..''••.• " 2~) (''-'• I,;*'" Tj) • J i,r,i,3+l ('V,; ''<■•' ~ TJ ) ( r .'.., • : r ',/ • ' ~ :j ) ) ( 22 ) 



i=0 3=0 

arid thr» int.prartinn with an PYtp.rna.l fiplH 



oo oc 



+ 1 



(23) 



tf/=££M<A^)=££>> 

i=0 j=0 V ' i=0 j'=0 

Formulas (|l9|), (p2|), (|23| ) realize the fermionic representation of the spin-i isotropic Heisenberg model on a 
square lattice pi), (0). 



IV. THE JORDAN— WIGNER TRANSFORMATION IN TWO DIMENSIONS (Y. R. WANG, 1991) 

Let us turn back to the Jordan- Wigner transformation in one dimension (|b|) using this case as a guideline and 
define a particle-annihilation operator as 



ck = e lQi S; , o-i = ^2 B ii n i> 

where By is the c-number matrix element, n- } = d^d$. A particle-creation operator is given by 

df =sfe iai =e iQi s+, 
whereas the inverse to Eqs. (p4j), ( p5| ) formulas read 

s7 =c iai d i = d i e [ai , 



(24) 



(25) 



(26) 



st =c- lai dt =dte~ lai . 



(27) 



We want the introduced operators d + , d to obey the Fermi type commutation rules. They are indeed the Fermi 
operators at the same site due to Eq. (||). Consider further two different sites i ^ j. Assuming that d + , d are 
Fermi operators one finds 



Since 



e -lnjBy d + e ln,B Jldj _ d . e -ln jBii d + 



(28) 



e~ in > B * = 1 + (e~ iB » - l) nj, e iniB - = 1 + (c is - - f ) m 
one can easily proceed in calculation of the r.h.s. of Eq. (^8|) finally arriving at 

= (f + c- [B -e [B Ad+dy 



(29) 



The result equals to (as it should for the commutator of the operators s + and s attached to different sites) if 



(30) 



Thus we assume that By 's in ( |24| ) , (|25|) , ([26|) , ( |27| ) satisfy relation (|30|) that yields the spin commutation relations 
for s + , s (H), @ if d + , d are Fermi operators. 

Y. R. Wang suggested! the following choice for Consider two complex numbers 



(31) 



and 



Tj = jx+ijy, 



(32) 



which correspond to the sites i = i x n x + i y n y and j = j x ii x + j y n y , respectively. Here n x and n y are the unit 
vectors directed along x and y axes, respectively. Assume that 

Bu = arg (73 - n) . (33) 

Evidently 

c iBji _ e iarg(ri-Tj) = e i(axg(T 3 -n)±n) _ _ e iB i3 ^ ^ 

that is the required condition (|3(i|). Since ri — 73 = \t{ — 73 |e iar § < -' ri ~ T - i ' Eq. ( [33]) can be rewritten in the form 

By =1111111(13-71), (35) 
and hence the introduced transformations (p4j), (|25|), (|2"6|), (|27j) contains 



«i = Imln (rj — T;) (36) 

j(#i) 



It is worth noting that the transformation of M. Azzouz ( |14| ) can also be written as Eqs. (|26j), ( pTj ) with 

By = it (6(4 - j,) (1 - + , (-)</, - j y ) (1 - <^J) (37) 



where Q(x) is the step function, and a necessary condition for having spin to fermion mapping (30) fulfilled. 

The advantage of the choice of Y. R. Wang (|36| ) becomes clear when one tries to introduce an approximate 
treatment of the transformed Hamiltonian. After inserting (|2^), ( |27| ) into ([!]), (|^) one finds 



+Jy (V* - J) ( d+di -l))+-£h (d+di ~) . (38) 



Note that 



Q!j — a\ = ^ dr • A(r) 



(39) 



where A(r) = V r a r and hence 



A(r) = J^T [ Yl "r'lmln (r' x -r x +i(r' y - r y )) j n 2 



+ I Z ^Tmln (r' x - r x + i (r' y - r y )) ] n B 

V \r'(#r) 



r'(#r) 

The crucial approximation to proceed is to make the change in (Boh 



(r'-rr 



n r - (n r ) = (4) + 1^1 (41) 



Here the angular brackets denote the thermodynamical canonical average with the Hamiltonian (0). Thus, we have 
postulated the mean-field description of the phase factors in (|3^). A similar treatment was adopted in Refs. 17 ll| 
Apparently, assuming further in that (sj) =0 one should suppose that h = (although in Refs. [n],[n] the 
uniform magnetic field was included into the Hamiltonian) . In principle Eq. ([ll]) simplifies the problem drastically 
since one faces a tight-binding spinless fermions on a square lattice. However, in practice it is hard to proceed 
because of nonuniformity of that model. 

The Hamiltonian (^8|), (|39|), j40|), (f4l| ) describes the (charged) spinless fermions moving in a plane in an external 
uniform (classical) magnetic field which is perpendicular to the plane. Due to ( pl[ ) (rv) can be taken out from the 
summation and in the continuum limit the vector potential of the field A(r) can be written0 in the following form 

^ n z x (r' - r) 

r'(/r) 



A(r) = -{u r ., ^ ; j- 

(r — rj 



'^n z xr, (42) 



where So is the area of the elementary plaquette in the plane (see Fig. 5) 




FIG. 5. Fermions in the magnetic field, which appears within the mean-field treatment of the phase factors in the 
Hamiltonian (|3^). 

and L x = L y = L — > oo. The corresponding magnetic field H(r) immediately follows from Eq. (|42|) 



H(r)=rotA(r) = (n r ,> — 

■30 



JL JL JL 

dx dy dz 

-r v r x 



2tt 

(n r ,)—n z , 



(43) 



so that A(r) = |H(r) x r. The flux per elementary plaquette equals to 



*n = TTM • ,<?nn_ = 9.Tr/n.,\ = tt 



The vector potential A(r) $4 ) shown schematically in Fig. 5 is not convenient to be concerned with. Because 
of the gauge invariance one may perform a gauge transformation introducing a new vector potential A(r) which 
yields the same flux per elementary plaquette $o = Namely, assume that A(r) is such that 

ri+i,j 

a i+ i d - Oij = dr • A(r) = tt, (45) 

J i,j 

whereas 

O-i+lj+l - ati+l,j = Oii,j+l - "i+lj+l = Oii,j ~ a »j'+i = (46) 



(see Fig. 6). From Eqs. (J45J) , (|4J) one finds that 

j> dr • A(r) = tt, (47) 

and on the other hand 

' dr • A(r) = / dS ■ rotA(r) = / dS • H(r) = <J> , (48) 



and hence the flux per elementary plaquette $o remains equal to tt. 

Let us turn back to the Hamiltonian (|38|). Only now we are in position to proceed with the statistical mechanics 
analysis. Within the frames of the introduced mean-field treatment of the phase factors (f45|), ( pi| ) the Hamiltonian 
( |38|) can be rewritten as 

H = E \ Ji -j ( d t d i - *4) + H *+ H f> ( 49 ) 

(iJ> 

where even for initially uniform lattice J = J' = Jj_ = J' ± in the Hamiltonian Hxy one has 

Ji,j;i+l,j = — J> Ji,j;i,j + 1 = J: Ji+l,j\i+2,j = J, Ji+l,j;i+l,j+l = J (50) 

etc. (see Fig. 6). 
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As can be seen from Eqs. (49), ( p0|) the isotropic XY model can be examined now without making any other 
additional approximations since it corresponds to a model of tight-binding spinlcss fermions on a bipartite square 
lattice (see Section VI). Conversely, the Heisenberg model requires further approximations to proceed because of 
the interaction between spinless fermions (see Section VII) . 



V. THE J ORD AN— WIGNER TRANSFORMATION IN TWO DIMENSIONS (E. FRADKIN, 1989) 

The Fermi-Bose correspondence in two dimensions, i.e., the Jordan- Wigner transformation for two-dimensional 
spin-i systems on a lattice, was discussed even earlier!, however, without applications to the theory of concrete 
spin models. Consider a system of spinless fermions, i.e., a matter field, o(r) on the sites of a square lattice (Fig. 
7) 




FIG. 7. Towards the Hamiltonian (|l|). 
and gauge field A x (r), A y (r) on the links of the lattice (Fig. 7). The Hamiltonian of the system is 

H = tJ2 (<x + (r)e U * (r) a(r + n x ) + a + (r)e iA « (r) a(r + n y ) + h.c.) (51) 

r 

with the constraint 



a+(r)a(r) = 9 (A y (r + n x ) - A y (r) - A x (r + n y ) + A x (r)) 



(52) 



where 9 is the parameter which will be defined later. E. Fradkin showecfl that the gauge field can be eliminated at 
the expense of a change in the commutation relations of the matter field. Namely, introducing^ the Jordan- Wigner 
operators S(r), a + (r) which obey 



-(r')a + (r) = 5 r ,, r - e^a + (r)&(r') 



l _ 



with is = 7r the Hamiltonian ( |5l| ) is 



H = t^2 (a+(r)a(r + n x ) + a+(r)a(r + n^) + h.. 



(53) 



(54) 



Eq. ( p4|) is easily recognized as the Hamiltonian of the spin-i isotropic XY model on a square lattice with exchange 
interaction between nearest sites J = 2t and the correspondence s+ = a + (r), s~ = 5(r), = 5 + (r)a(r) — \. 

Further discussions on the extension of the Jordan- Wigner transformation to three- or more-dimensional cases 
can be found in Ref. 



VI. 2DS=i ISOTROPIC XY MODEL 



_• 



Let us show how the two-dimensional Jordan- Wigner transformation with the mean-field treatment of the phase 
factors can be used in the theory of spin models. As a result we come to an approximate approach to the study of 
2D quantum spin models. 

We begin with the 2D isotropic XY model (Hz — Ht — 0) considering for concreteness the case of J = J', 



J± = J' ± . In accordance with (p9[), (50) we start from 



1 v x 

HxY = 2 £ ^ ( J ( d ti di +^ - d ^ d t+i,j) + J± {di tj d id+1 - d itj df j+1 )) 

i=0 j=0 

= 2 J (• ■ ' - a ti & i+i.i + a iJ b t+i,j + b t+i,j a i+2,j - b i+ ijaf +2tj + ■ ■ •) 
+\j± {■■■ + af tj bi !j+ i - ai,jbf tj+1 + b+ +lj a l+1 . j+1 - b i+ld af +ld+1 + ...), (55) 

where, to emphasize a bipartite character of the square lattice that appeared, we have introduced notations o<j = 
di,j, l>, . | , = di+ij etc.. 

Notwithstanding the fact that we have an approximate theory, it contains the exact result in ID limit if one puts 
in Eq. (^5|) either J± — or J = coming to a system of noninteracting chains extended in either horizontal or 
vertical direction, respectively. In the latter case Eq. (|5^) corresponds to a system of noninteracting chains each 
with the XY part of the Hamiltonian (rU 



oo ^ 

Hxy(i) = o J ^ RtAj+i - d i,i d i,j+\) ■ ( 56 ) 



2 

j=0 



In the former case one gets a system of noninteracting chains each with the Hamiltonian 

HxyU) = £ 2 J WA+hi - dijdi. id ) (57) 

i=0 

and to recover the ID limit explicitly the transformation 

<, = e im/, */+, d i+hj = e-^Vi+i, • • • , ^o = 0, Vi+i = ^ + i (58) 

(e.g., d+ 3 = /+, d+ d = -/+, d^j = -/ 2 +, d£ d = / 3 + etc.) should be performed. 
After the Fourier transformation 

di , = ; 1 V e^+^dk k (59) 

etc. or in short 

di = ^Vc iki d k , d+ = ^Ve- iki rf+ i^Ve^i J= Vc iki rf+ 

Va ^ ^ k Va ; Va 7 ; 



2tt N x N x N x 2tt N y N y N v 

n x , n,- — , — + !,..., — -!, n y , n y -- — ,-— + !,..., — -!, 

/a. r/+ "!• = . Jd,. . rf,.A = . d+ "!• = n rfin^i 



(N x , N y are even), the Hamiltonian (|5q) becomes 



H X y = ^ X! ( iJ smfc:E ( & k ak - a k fo k) + J± cos k y (b£b k - a^a k )) 

k 

= o ^ l-^kl (cos 7 k (6+6 k - a+a k ) + isin7 k (&k a k - a^6 k )) 



where we have used the notations 



b k- d k^k y ' °k - d k:c ±7T,k v ±7T, 

-Ek = </l cosfc a + Usmk x = |£' k |e l7k , 



Ek | = y J\ cos 2 k y + J 2 sin fc^ , cos 7 k 



J± cos ky 



sm7 k 



J sin k x 



The Hamiltonian (61) can be rewritten as 

H X y =2J'|-Ek| (cos7 k (fej^k - a k a k ) + isin7 k (b£a k - a£b k )) 

k 

where the prime denotes that k varies in the region shown in Fig. 8b. 

k y a 



K k 



(61) 



(62) 



(63) 




K k 



FIG. 8. The region in which k varies in the sum in Eq. ( |si| ) (see Eq. (|6C|)) (a) and in Eq. ( |63| ) (b). 
At last one introduces the operators 

7 k 7 k 7 k 7 k 

a k = cos —b k + isin -yak, /3 k = sin —b k - icos y ak, 

{"ki,"^} = <5ki,k 2 , {/3ki,^k 2 } = ^ki,k 2 , {a kl ,#£J = {Aci," k2 } = 0, 



(64) 



etc. to get from Eq. (|6^) the final form of the 2D spin-^ isotropic XY model Hamiltonian in fermionic language 
within approximation (fill ) 



H XY = ' A k K« k - , 



A k = I S k | = \l J 2 sin 2 fc^ + J\ cos 2 fc. u > 



(65) 



It is easy now to calculate the thermodynamic functions of the spin model which correspond to Eq. (pa). For 



p-sramnlp tlnp Q'rmmrl sta.tp pnprcrv npr site is 



1 r r /j2 sin 2 fc + j2 cos2fe 

2 j_ v 2tt 2^ v -l » 



(66) 



In ID limit ( J ± = or J = 0) Eq. (pFJ) becomes 



at = -2 /_ ^ jsmkxi = ~2 y_ ^ cos ^' = - V (67) 

that is the well known exact result. 

It is remarkable to note that the Hamiltonian ( p5| ) (and hence the ground state energy per site (|66|)) arises also 
for the 2D spin-i Heisenberg model and is known as the uniform flux ((n r ) =4) solution!. 



VII. 2DS=i HEISENBERG MODEL 



2 



In the remainder of the paper we consider the 2D spin-i isotropic Heisenberg model ([!]), (0) treating the phase 
factors which appear after making use of the 2D Jordan- Wigner transformation within the frames of the mean-field 
approximation ([fl]) . The Heisenberg model besides the Hxy term (^H|) includes the Hz term which contains 

d+ckd+di - \d+di - ±dfd s + \. (68) 

Thus, the Jordan-Wigner spinless fermions interact and further approximations are required. The first term in 
can be changed by 

d+did+di -> d+di(dfdj) + (d^d^dfdj - (d+diXdfdj) 
+d+dj{didt) + (d^dj)didf - (dfd^ididt) 

= d tM n j) + d t d i( n i) - ( n i)( n i) 
+d 1 +d j A iJ e ie - - didfAj-J^ •> + AyA J ,ie i ^J + °i- ,) ) (69) 

where we have introduced the notation 

(didf) = A u e^ = -(dtdi). (70) 
In accordance with (^9|) there may be four ways to treat the Ising interaction, i.e. assuming either 

Ay=0, m = ( S |) = (nj)-i = (71) 

or 

Aij ^ 0, m = (72) 



or 



or 



Ay = 0, m ^ (73) 



Ay^O, m^O. (74) 
The first possibility ([n]) yielding the uniform flux solution was considered in the previous Section (Eqs. fl65|), 



Let us consider the second possibility (|72|). In such a case the Ising term becomes 

H Z = £ J,j - -\j.ic'"' 1 </;<// + A iJ A j il - i "^ i "i: ) , 

where to get the in-phase flux solution one puts 

glfli.ijl+l.j _ gW*+l,i;»,j _ e 19 i+l,3;i + 2,j _ glSi+2,3i»+l,3 _ J 

Aij;ij + 1 = AiJ + l;iJ = Ai+l,j;«+l,j + l = Aj + lJ + + J = P, 

e W *.ii*.i+l = gW«,i+lii,j — gl*i+l.i!i+l,J+l — g ie i + i.j + i;i + i,i _ 2 

(see Fig. 9) 

















G 












j+ij+i 
p 










i+hj 















FIG. 9. Towards the in-phase flux solution for the Heisenberg model; the values of Aye id are attached to the bonds. 

and the parameters Q and P are calculated self-consistently (see Eqs. (j30|), ( |8l| ) below). Now the Heisenberg 
model Hamiltonian only slightly differs from that of the XY model (^35|) becoming 

H = - J (1 + 2Q) (. . . - a+-6 l+1 j + a itj bf +ltj + bf +l j a i+2j - 6 i+lj -a+ 2J + . . .) 

+- J ± (1 + 2P) (. . . + af j b i;j+1 - a itj bf j+1 + b+ +l j a i+1 ^ +l - b i+hj af + + . . .) 

+7VJQ 2 + NJ ± P 2 . (77) 

Acting along the line described in the Section VI one hnds that 

H = J2 ' A k (<*k«k - p£0y) + NJQ 2 + NJ ± P 2 , 

k 

A k = \jj 2 (1 + 2Q) 2 sin 2 fc^ + Jl(l + 2P) 2 cos 2 /c y > 0. (78) 
The ground state energy per site is given by 

f = ~\ fjt f v (1 + 2Q)2 Sin2 K + Jl(1 + 2Pf C ° S2 *» 

+J0 2 4- J i P 2 r7Ql 



(75) 



(76) 



n i r dk x r dk y j S in 2 Mi + 2Q) 



where the parameters Q and P are determined from the conditions -§q^t — and -§pj$- = 0, i.e., 

1 + 

1 J 2 (1 + 2Qf sin 2 fc x + J 2 (1 + 2P)^ cos 2 k 

and 

1 dfc^ [* dk y Jj_cos 2 k y (l + 2P) 



2 -' " " >rr 7 ' 2; ' 7 J 2 (1 + 2Q) 2 sin 2 k x + J 2 (1 + 2P) 2 cos 2 fc y 



Treating the Ising term (|69| ) within assumptions (|73|) or (|74| ) one assumes the Neel order with the sublatticc 
magnetizations m and — m 



(m,j) = (n l+1J+1 ) = . . . = m + i, (n^+i) = (n<+i, 3 -) = . . . = -m + ~ (82) 



(see Fig. 10). 
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FIG. 10. The Neel order for the Heisenberg model; the values of magnetization (sf ) are attached to the sites. 

(The given assumption, by the way, apparently contradicts the mean-field approximation ( pl| ) for the phase 
factors; this inconsistency, however, to our best knowledge has not yet been discussed.) Therefore, case ([73]), i.e., 
the uniform flux with the Neel order, is called the Neel flux solution and case (f74|), i.e., the in-flux with the Neel 
order, is called the in-phase Neel flux solution. In the latter case the Hamiltonian reads 

H = -J (1 + 2Q) (. . . - a+-&. 1+ i j + ciijbf^j + bf +lj a i+2 ,j - b i+ld af +2J +...) 

+-J± (1 + 2P) (. . . + af j b itj+1 - a^b+, J+l + b+ +hj a i+ i, j+ i - b l+1 . ja + +lj+1 + . . .) 

+J (. . . - mafjOij + mbf +ld b i+1 j + ////-,. , ,. , - ma+ +2j a. 1+2j + . . .) 
+ Ji (• . . - + mi)f j+1 b iij+ i + mbf +ld b i+ i tj - ma+ +lj+1 a i+hj+1 + . . .) 

+NJQ 2 + NJ ± P 2 + N (J + J x )m 2 . (83) 

Performing the Fourier transformation ( |60| ) one gets 

H = ^Y1 ((2 (J + J±) m + i J (1 + 2Q)smk x )b+a u 



+ (2 ( J + J x ) m - i J (1 + 2Q) sin a^6 k 
+ Jx (1 + 2P) cos fc y (fe+6 k - a+a k )) 
+NJQ 2 + NJ^P 2 + N(J + Jj_) m 2 
^ ' ((2 (J + Jj_) m + iJ (1 + 2Q) sin fc^) 6+a k 

k 

+ (2 ( J + Jj_) m - iJ (1 + 2Q) sin k x ) a+6 k 
+ Jx (1 + 2P) cos fc y (6+fe k - a+a k )) 

+7V JQ 2 + AT J ± P 2 + AT ( J + J±)m 2 , (84) 



where b£ = k , a k = dk^±-K,k y ±ni etc.. Introducing the operators 



2 ( J + J ± ) m ± i J (1 + 2Q) sin fc^ = \Ja{J + J ± f m 2 + J 2 (1 + 2Q) 2 sin 2 k x e ±i<5k (85) 
and then the operators 

w k r . w k _ Wkr ^k ~ 

«k = cos — 6 k + sin -y a k , /j k = sin — b k - cos -y «k; 



Jx (1 + 2P) cos/c„ 
cosa; k = — — -, sma; k — 



^4 (J + Jx) 2 m 2 + J 2 (1 + 2Q) 2 sin 2 k x 



' |Pk| 



|£ k | = y 4 ( J + Jx) rn 2 + J 2 (1 + 2Q) sin 2 k x + J 2 (1 + 2P) cos 2 k y (86) 
one gets the final form of the Heisenberg Hamiltonian in fermionic language 

H = ' A k K«k - /3 k /3k) + A^JQ 2 + JV JxP 2 + A r ( J + Jx) m 2 (87) 

k 

with A k = \Ek\ > defined by Eq. @. 

The ground state energy per site which follows from ( |S7| ) reads 

lf = -\f_l£f_ ( J + Jx) 2 ™ 2 + J 2 (1 + 2Q) 2 sin 2 + J 2 (1 + 2P) 2 cos 2 k y 

+JQ 2 + J ± P 2 + (J + Jx) m 2 , (88) 

where the values of the introduced parameters are determined by minimizing ^9- (pq) with respect to Q, P and m 



N 

f , Jsin 2 fc ,(l + 2Q) (8Q) 



P= W " , ; (90) 



(27T) 2 , 


/ 4 (J- 


- Jx) 2 TO 2 " 


h J 2 (1 + 2Q) 2 sin 2 k x 4 


-J 2 (l- 


h 2P) 2 cos 2 fc y 


dk 






Jx cos 2 fcj, (1 + 2P) 






(27T) 2 , 


/4(JH 


- Jx) 2 fn 2 - 


1- J 2 (1 + 2g) 2 sin 2 fc x 4 


-J 2 (l- 


h 2P) 2 cos 2 fc y 


dk 






2 (J 4- Jx)m 






(27T) 2 y 


A (J 4 


- Jx) 2 m 2 H 


- J 2 (1 4- 2Q) 2 sin 2 k x 4 


•Ji(n 


- 2P) 2 cos 2 k y 



2m = / ? v (91) 

7 ( 2?r ) ^/4(J + Jx) 2 m 2 + J 2 (1 4- 2Q) 2 sin 2 fc x 4- J ± (1 4- 2P) 2 cos 2 fc,, 

The Neel flux solution (sec Eq. @) is given by Eqs. @, @, @, (§§), @ in which P and Q are equal to 



zero. 



VIII. THE 2DS=i HEISENBERG MODEL: A COMPARISON OF SOME RESULTS 



In the final Section we want to list out some problems of the two-dimensional spin models theory which were 
attacked with the exploiting of the 2D Jordan- Wigner transformation. 

The 2D antiferromagnetic Heisenberg model with J = J 1 = J± = J' ± (without external field) was considered 
m Refs. gMg. The main results obtained concern the ground state energyi'B, the specific heatfl and the Raman 
spectrumQu. A comparison with some experimental data for La2Cu04 was given. 

In Ref. |§| the effects of the interchain interaction on the one-dimensional spin-^ antiferromagnetic Heisenberg 
model were examined. For this purpose the Heisenberg model with J = J' , J± = was considered and the 
in-phase Neel flux solution was analysed at zero temperature. The author found that the one-dimensional limit 
is singular, i.e. the staggered magnetization m/0 (2m = 0.513) when J± — > +0 (although we know from exact 
results that for the antiferromagnetic chain m = 0). In the other limiting case J = J± the theory based on the 
fermionization procedure yields 2m = 0.778 (the spin wave result is 2m = 0.6; more accurate calculations predict 
m = 0.3074 (see Ref. |l|)). The result of Ref. || stays somewhat separately in the estimate of the value of ^j- at 
which the staggered magnetization appears. Different theories predict ^j- from to 0.2 (for details see Ref. pj| ). 
Ref. H predicts ^j- = 0, moreover 2m jumps from zero to 0.513 for any infinitesimally small 

The antiferromagnetic Heisenberg model on a ladder within the frames of the in-phase Neel phase solution 
(the ground state energy, the singlet-triplet energy gap) was discussed in Ref. |t| The effects of the interladder 
interaction and magnetic field on the susceptibility at nonzero temperatures were studied in Ref. [ll]. 

A consideration of the spin-Peierls state under magnetic field was reported in Ref. [l0|. A study of the stepped 
spin-Peierls transition for the quasi-one-dimensional XY and Heisenberg models using the 2D Jordan- Wigner 
transformation was reported in Refs. |l2|,[l3| The 2D Jordan- Wigner transformation was applied for a study of the 
zero temperature spin-Peierls transition for the quasi-one-dimensional XY and Heisenberg systems in Ref. [l4|. 
In particular, the phase diagram between the dimerised and uniform states in the parameter space of interchain 
interaction and spin-lattice coupling was constructed. 

Many more problems may be considered within the frames of the 2D Jordan- Wigner fermionization approach. 
Probably, the 2D Jordan-Wigner transformation should be of more use for the 2D spin-i XY models since for 
such models no further approximations (except the mean-field-like treatment of the phase factors) are required. 
Besides, a more sophisticated treatment of the phase factors is desirable. 

The author acknowledges the kind hospitality of the Max Planck Institute for the Physics of Complex Systems, 
Dresden in the end of 1999 when the paper was launched. He thanks J. Richter (Magdeburg University) for the 
kind hospitality in the summer of 2000 when the main part of the paper was prepared and discussed. Special thanks 
go to T. Krokhmalskii, T. Vcrkholyak, O. Zaburannyi and V. Derzhko for helpful comments and suggestions. 
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